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The AGS equations are solved for n 3 H and p 3 He scattering including the Coulomb 
interaction. Comparison with previous work confirms the accuracy of the calculation and 
helps clarify a number of issues related to the n 3 H total cross section at the peak of the 



resonance region, as well as an A y deficiency in p 3 He. Calculations are fully converged in 
terms of NN partial waves and involve no uncontrolled approximations. 



1. INTRODUCTION 

VO ' 

In spite of the progress that has taken place over the years on the solution of the three- 
nucleon scattering problem in terms of ab initio calculations, there are still some open 
problems that defy our understanding. They are the A y puzzle in elastic Nd scattering 
and the space star anomaly in Nd breakup, just to name a few. These problems persist 
even when the Coulomb interaction is added [Tf2"f3"] and seem to be insensitive to the 
choice of realistic 2N + 3N force model. As it has already been highlighted in the past [U 
EEL the four-nucleon scattering problem reveals further discrepancies between theory and 
experiment that need clarification in terms of improved calculations using modern force 
models and efficient numerical algorithms that allow for a numerically converged solution 
of the Alt, Grassberger and Sandhas (AGS) equation [7] for the transition operators. 

In four-nucleon scattering the Coulomb interaction is paramount to treat p 3 He, to 
separate the n 3 He threshold from p 3 H and at the same time avoid a second excited state 
of the alpha particle a few keV bellow the lowest scattering threshold [SJ. Given the 
success recently achieved in including the Coulomb interaction in pd elastic scattering 
and breakup [2|3] . we generalize the proposed method to the four-nucleon system and 
study p 3 B.e below three-body breakup threshold. 

The work we present here constitutes the first effort to design and construct a new 
scattering code that may bring the treatment of the four-nucleon system to the same 
degree of accuracy and sophistication as we already have in the three-nucleon system. 
The symmetrized AGS equations for four identical nucleons are given by 

U 11 = -(G t Go)- 1 P 34 - PuU G t G U 11 + UG tG U 2 \ (la) 
U 21 = (Go t Go)' 1 (1 - P34) + (1 - Pm)U Go t G U n , (lb) 

where U 11 (U 21 ) is the transition operator for 1+3^1+3 (1+3^2+2), t the NN t-matrix 
and Go the four free particle Green's function. The operators U and U are respectively 
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1 + (3) and (2) + (2) subsystem transition operators given by 

U = PGq 1 + PtG U, (2a) 
U = PG Q l + PtG U, (2b) 

where the permutation operators P and P are given by P = P12P23 + P13P23, and P = 
-Pi3-f24- Defining and | ^2) the initial/final (1 + 3) and (2 + 2) states that satisfy the 
equations 

\<h) = G tP\<i> 1 ), (3a) 
\<p 2 ) = GotP\<p 2 ), (3b) 

the matrix elements {(p a \T a0 \(j)p) = S a ((p a \U a ^\(f)f3) with Si = 3 and S 2 = 2 lead to the 
transition amplitudes from which one calculates observables. After partial wave decom- 
position we solve a three-variable integral equation with no approximations beyond the 
usual discretization of continuous variables in a finite mesh. 



2. RESULTS 

In order to calibrate the accuracy of our new numerical algorithm we calculate first 
the negative eigenvalues of the kernel that correspond to the 4 He binding energy. The 
Coulomb interaction is included through screening and charge dependence in the NN 
interaction taken in isospin T = alone. In Table 1 we show the convergence of the 
results for the binding energy in terms of two-nucleon total angular momentum I. The 
calculation with I < 6 includes 1057 channels and clearly shows the accuracy of the 
numerical method we use. Results of ref. [12] for CD Bonn include total isospin T = 1 
and T = 2 states, whereas they are neglected in our calculations. 



Table 1 

a-particle binding energy as function of the total two-nucleon angular momentum /. 
Predictions for AV8' do not include Coulomb. 





/ < 1 


I < 2 


I< 3 


I < 4 


I< 5 


I < 6 


refs. pjnnnni2] 


AV8' 


23.08 


25.16 


25.69 


25.85 


25.90 


25.91 


25.90 - 25.94 


AV18 


22.30 


23.75 


24.15 


24.20 


24.23 


24.24 


24.22 - 24.25 


CD Bonn 


25.03 


25.95 


26.07 


26.10 


26.11 


26.11 


26.13 - 26.16 


INOY04 


28.68 


29.09 


29.10 


29.11 


29.11 


29.11 


29.11 



Next we show results for n 3 H scattering by solving Eqs. (flaj) and (llbl) in isospin T = 1. 
The phase shifts at E n = 3.5 MeV are shown in Table [2] together with the total cross 
section at. The present work confirms previous calculations by the Grenoble and Pisa 
groups (second and third lines respectively) compiled in ref. [6] together with the results 
by one of us (fourth line) where the AV18 interaction was represented in rank one. The 
small differences between the phases in the first three lines are probably due to higher iViV 
partial waves (J < 4) that are included in the present calculation together with £ y , £ z < 4. 
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Table 2 



n 3 H phase shifts and mixing parameter at E n = 3.5 MeV together with total cross section. 





0+ 


o- 


1+ 


1- 


2~ 








(%) 


( 3 Po) 


( 3 Si) 


( 3 Pi) 




5 


( 3 P 2 ) 






AV 18 


-66.07 


20.72 


-58.44 


40.08 


20.71 


-44.50 


42.48 


2.329 


this work 




-66.5 


20.9 


-58.5 


37.3 


20.7 


-43.5 


41.0 


2.24 


ref. [6] 




-66.3 


20.6 


-58.7 


38.6 


20.5 


-45.5 


40.1 


2.24 


ref. [6] 




-63.7 


27.5 


-58.3 


44.8 


24.7 




44.2 


2.51 


ref. [6] 


CD-Bonn 


-64.63 


18.97 


-57.41 


39.44 


20.19 


-44.95 


42.45 


2.282 


this work 



In spite of the improved accuracy, both AV18 and CD Bonn interactions lead to cross 
sections in the resonance region that fail to reproduce the experimental data. This is 
shown in Fig. [TJ 




Figure 1. Total cross section for n 3 H scattering versus center of mass energy for different 
realistic interactions. 



Finally we use the method developed for pd elastic scattering [2|3] to include the 
Coulomb interaction between the protons. This means using the two-potential formula 

(4>M Tll \M) = (PfMPi) + Km {^(P/K^P/I [Tg, - t R ] |&a)2£'(ft)} , (4) 

where the first term is the long range Coulomb amplitude between the proton and the 

center of mass of 3 He and the second term represents the Coulomb modified nuclear short 

range amplitude that results from the difference between the matrix elements of the AGS 

operator calculated with screened Coulomb between the three protons and the two 

_ i 

body p 3 He screened Coulomb t-matrix tn after renormalization with Z R 2 . 
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Figure 2. Convergence of p 3 He observables at E p = 4 MeV with the screening radius R 
of the Coulomb potential. 




Figure 3. Differential cross section and nucleon analyzing power for p 3 He scattering at 
E p = 5.54 MeV using CD Bonn (solid) and AV18 (dotted) potentials plus the Coulomb 
interaction between all three protons. The dashed line corresponds to CD Bonn alone. 
Experimental data are from refs. [l~3~yT4"] . 



The convergence of the calculation with the screening radius R is shown in Fig. [2j For 
R = 12 fm we get converged results for both observables. Further results are shown 
in Figs. [3] and H] for the da/dfl, A y , C xx and C yy at E p = 5.54 MeV. The calculations 
show that there are large Coulomb effects, and greater sensitivity to the NN interaction 
compared to what is observed in pd scattering at low energy. The calculated analyzing 
power A y shows a large discrepancy with the data but C xx and C yy are in reasonable 
agreement with experiment. 
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Figure 4. Spin correlation coefficients for p 3 He scattering at E p = 5.54 MeV. Curves as 
in Fig. [3j Experimental data are from ref. |14j . 
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